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We use the Newtonian limit of a general scalar-tensor theory around a background field to study 
astrophysical effects. The gravitational theory modifies the standard Newtonian potential by adding 
a Yukawa term to it, which is quantified by two theoretical parameters: A, the lenghtscale of the 
gravitational interaction and its strength, a. Within this formalism we firstly present a numerical 
study on the formation of bars in isolated galaxies. We have found for positive a that the modified 
gravity destabilizes the galactic discs and leads to rapid bar formation in isolated galaxies. Values 
of A in the range « 8 - 14 kpc produce strongest bars in isolated models. Then, we extent this work 
by considering tidal effects due to interacting galaxies. We send two spirals to collide and study 
the bar properties of the remnant. We characterize the bar kinematical properties in terms of our 
parameters (A, a). 
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Usamos el h'mite newtoniano de una teon'a escalar-tensorial general alrededor de un campo del 
fondo para estudiar efectos astroHsicos. La teon'a gravitacional modifica el potential newtoniano 
estandar, agregandole un termino de Yukawa, el cual se cuantifica por dos parametros teoricos: A, 
la escala de longitud de la interaccion gravitacional y a, su intensidad. Dentro de este formalismo 
primero presentamos un estudio numerico de la formacion de barras en galaxias aisladas. Encon- 
tramos que para a positiva la gravedad modificada desestabiliza a los discos galacticos y Ueva a una 
rapida formacion de la barra en galaxias aisladas. Valores de A en el rango de « 8 - 14 kpc producen 
barras mas pronunciadas en los modelos de galaxias aisladas. Despues, extendemos este trabajo al 
considerar efectos de marea debido a galaxias que interaccionan. Enviamos dos galaxias espirales a 
colisionar para estudiar las caracten'sticas de la barra del remanente formado. Caracterizamos las 
propiedades cinematicas de la barra en terminos de nuestros parametros (A, a). 

Descriptores: Formacion de la barra; interaccion de galaxias; teorfa escalar-tensorial 
PACS numbers; 04.50.-|-h, 04.25.Nx, 98.10.+Z, 98.62.Gq, 98.62.Js 



I. INTRODUCTION 

In recent years there have been some attempts to ex- 
plain observed gravitational effects that imply the exis- 
tence of dark matter (DM) and dark energy (DE) . These 
effects go from galactic to cosmological scales and we still 
do not know the nature of DM and DE, or if they are re- 
lated to each other, or even if DM is unique or there is 
a set of dark multi-components. One possibility is that 
scalar fields (SF) play a role in the modified dynamics. 
Particularly, SF are often applied in cosmology to ac- 
complish the Universe's accelerated expansion, inferred 
from Supernovae la redshifts, CMB Doppler peaks mea- 
surements, large-scale structure surveys and cosmological 
simulations 1, 2, 3, 4, 3^^,0,0- For a review of all these 
topics see for instance [9|. 

The way in which SF couple to gravity is also unknown, 
simply because there is a lack of a unique fundamental 
theory that explains the intricate relationship of matter 
and its gravitational background. One possibility is that 
SF are coupled non-minimally to gravity, at a Lagrangian 
footing, as it happens when string theories are compact- 
ified to four space-time dimensions fl^. The resulting 



effective theory is a scalar-tensor theory (STT) of grav- 
itation, that can generically be described by arbitrary 
scalar functions, apart from the geometrical part 
In the past we have studied different effects of this type 
of theories in cosmology [H, [H, [H [H, [13, d, and 
more recently we have considered the Newtonian limit 
of STT and apply it to astrophysical phenomena. We 
have computed potential-density pairs for various halo 
density profiles [l^ and axisymmetric systems HO]. It 
was found that the augment of rotation curves and ve- 
locity dispersions depend on the parameters of the SF. 
On the other hand, in Ref. [2l[ we have computed the ef- 
fect of SF on the transfer of angular momentum between 
protogalaxies. In the present work, we pursue to study 
bars in spiral galaxies. Observations of spiral galaxies 
indicate that the presence of a bar is a common feature 
[2^ . Instabilities in isolated stellar and gaseous discs lead 
to bar formation; see [231 for pioneer studies and 0, [2^ 
for a modern view. The bar formation in isolated models 
has been widely studied both analytically and numeri- 
cally [e.g., m, [23, Sin, HO, [Ml,!!!, [Hi and it is studied 
using the above STT formalism in the first part of this 
paper. In the second part, we consider in dynamical ef- 
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fects of non-isolated systems which are found in clusters 
of galaxies. In this sense, it has been suggested that the 
observed bar in many spirals is the result of the gravita- 
tional interaction between two or more nearby galaxies. 
For instance, Nogushi has found that during the col- 
lision of two galaxies and between the first and the second 
closest approaches, the disc takes a transient bar shape. 
The gravitational interaction between the two galaxies 
gives rise to perturbations in the orbits of the stars that 
results in the formation of the bar. 

Bar formation in stellar discs depends upon various 
simultaneous effects. In the case of collisions, the 2D- 
simulations have shown that these factors are [3a |: rota- 
tion curve shape, disc-halo mass ratio, perturbation force 
and geometry. Additionally, simulations suffer from nu- 
merical effects such as low spatial and temporal resolu- 
tion, too few particles representing the system, and an 
approximate force model. These effects were studied by 
us m Rcfs. [25,^, where it was shown that specific 
parameter choices may change bar properties. Once nu- 
merical effects are controlled, we may investigate all the 
other model parameters, which in our case are (A, a). 

In the present paper we study the formation of bars 
as a product of both instabilities of isolated galaxies and 
as a result of the collision of two spirals in framework of 
STT. In particular, we consider a non-minimally coupled 
SF in the Newtonian limit (section |TT| and use the result- 
ing modified gravitation force in our SD-simulations. In 
this way, all collisionless particles mutually interact with 
the modified gravitational force. Then, we investigate 
isolated galaxies (section IIIip and head-on and off-axis 
impacts of two disc galaxies and the properties of tidally 
formed bars (section lIVp . We finally draw some conclu- 
sions (section IV)) . 



II. SCALAR-TENSOR THEORY AND ITS 
NEWTONIAN LIMIT 



A typical STT is given by the following Lagrangian 
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from which we get the gravity and SF equations. Here 
g^^ is the metric, CM{g^v) is the matter Lagrangian and 
Ci;((/)) and V{4>) are arbitrary functions of the SF. Thus 
the gravitational equation is 



7? -i„ /? - i 



SttT^. + ]^Vg^, + d,<j) (2) 



-Y^{^^,<^)f■g^,v + - g^J.v^i 

The SF part is described by the following equation 
4>V' - 2V 1 



[%T,T-u\d^f] , (3) 



where a prime (') denotes the derivative with respect to 
SF(0). 

In accordance with the Newtonian approximation, 
gravity and SF are weak. Then, we expect to have small 
deviations of the SF around the background field. As- 
suming also that the velocities of stars and DM particles 
are non-relativistic, we perform the expansion of the field 
equations around the background quantities (0) and 77^1^ . 
Even though the expansion of the above equations to first 
order is well known [s^, [H, [s^ , we explicitly show it in 
the appendix since our definition of the background field 
is (0) = G'^^(l-l-a), which is non-trivial, and this changes 
some constant terms in the equations. Accordingly, we 
obtain Eqs. ([AT6l) and (1X221) : 
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where p is matter density of DM or stars stemming from 
the energy-momentum tensor, Gn is the Newtonian grav- 
itational constant and a = 1/(3 + 2uj) is a constant, in 
which oj is the Brans-Dicke parameter fTl|, here defined 
in theories that include scalar potentials. Equations (j4|) 
and ^ represent the Newtonian limit of a set of STT 
with arbitrary potentials (y((/))) and functions uj{4>) that 
are Taylor expanded around some value. The resulting 
equations are thus distinguished by the constants a and 

TO. 

In the above expansion we have set the cosmological 
constant equal to zero since within galactic scales its in- 
fluence is negligible. This is because the average den- 
sity in a galaxy is much larger than a cosmological con- 
stant that is compatible with observations. Thus, we 
only consider the influence of luminous and dark matter. 
These matter components gravitate in accordance with 
the modified-Newtonian theory determined by Eqs. ^ 
and (O. The latter is a Klein-Gordon equation which 
contains an effective mass to term, whose Compton wave- 
length (A = h/mc) implies a length scale for the modified 
dynamics. We shall assume this scale to be of the order 
of tens of kilo-parsecs, which corresponds to a very small 
mass, to ~ 10~^^ eV. 

Note that Eq. ^ can be cast as a Poisson equation for 
^ = (l/2)(/ioo + ^/(0)), 



= AttGn p/(1 + a) 



(6) 



Thus, the modified Newtonian potential is now given by 



(7) 



Particular solutions, the so-called potential-density 
pairs [40|1, were recently found for the NFW's and 
Dehnen's density profiles and for axisymmetric sys- 
tems [2^. For point masses (of non-SF nature) the so- 
lution is well known [3^ |4l| and here is adapted to our 
definition of the background field, (0) — G^^(l + a): 
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where 

_ Gn \ - rus 
" - (i + „)Z.|r-r,|' 

with TTT-s being a source mass. The potential u is the 
Newtonian part and u\ is the SF modification which is of 
Yukawa type. The total gravitational force on a particle 
of mass rrii is 

F = -miV$Ar = TO,;a. (11) 

Thus, gravitating particles, that in our simulations are 
stars or DM particles, will feel the influence of Newtonian 
gravity (u) plus a SF force due to the term ux- 

The gravitational potential of a single particle arising 
from the above formalism is: 

<^>^^^^^{l + ae-^/^), (12) 
(1 + a)r 

For local scales, r ^ A, deviations from the Newtonian 
theory are exponentially suppressed, and for r ^ A the 
Newtonian constant diminishes (augments) to Gn /{^ + 
a) for positive (negative) a. This means that equation 
fulfills all local tests of the Newtonian dynamics, and 
it is only constrained by experiments or tests on scales 
larger than -or of the order of- A, which in our case is of 
the order of galactic scales. By contrast, if one defines 
< 4> >= ^/Gn, then the effective Newtonian constant is 
modified at scales r < A, and stringent, local constraints 
applies, demanding a to be less than 10"^" [4l|. This 
latter approach will not be considered here. 

Recently, the effect of STT has been investigated in 
different cosmological scenarios in which variations of the 
Newtonian constant are constrained from a phenomeno- 
logical point of view. For instance, [13] studied the influ- 
ence of varying Gat on the Doppler peaks of the CMBR, 
and concluded that their parameter (^ = G/Gn) can be 
in the interval 0.75 < ^ < 1.74 to be within the error bars 
of the CMBR measurements. In our notation this trans- 
lates into —0.43 < a < 0.33. However, this range for a 
has to be taken as a rough estimation, since these au- 
thors have only considered a variation of Gn, and not a 
full perturbation study within STT. The latter has been 
done by [11] , who found some allowed deviations from the 
Newtonian dynamics, that translated into our parameter 
is a = 0.04; however, a comparison with observations 
in not made. On the other hand, a structure formation 
analysis has been done in [3| , in which deviations of the 
matter power spectrum are studied by adding a Yukawa 
potential to the Newtonian. They found some allowed 
dynamics, that turns out to constrain our parameter to 
be within —1.0 < a < 0.5; but again a self-consistent 
perturbation study in general STT is missing. Thus, the 
above three estimates can be taken as order-of-magnitude 
constraints for our models. Note that even when it is not 



theoretically justified to take negative values for a, phe- 
nomenology admits them. In this work, however, we only 
consider positive values of a. 



III. ISOLATED GALAXY SIMULATIONS 

We use the standard procedure to construct a ga laxy 
model with a Newtonian potential described in [2^, |36[ . 
The galaxy consists of a disc, halo, and bulge and its 
initial condition is constructed using the Hernquist halo 
model (a Dehnen's family member with 7=1, see [3). 
To perform the SD-simulations we used the gbsph code 
(www. astro . inin.mx/mar/nagbody) modified to include 
the contribution of the scalar fields as given in the pre- 
ceding section. The forces were computed with a toler- 
ance parameter 6 = 0.75, and including the quadruple 
term. We use Barnes's model parameters and system 
of units [l^. The mass, length and time scales are set 
to 2.2 X 10" Mq = 1.40 kpc = 1 and 250 Myr= 1, re- 
spectively. In these units, the gravitational constant is 
Gn = 1. The discs scale height is zq — 0.007 and the half 
mass radius of the galaxy is located at R1/2 ~ 11 kpc. 

All isolated runs were performed with e = 0.015 {— 0.6 
kpc) and At = 1/128 for N = 40 960, and e = 0.008 
and At = 1/128 for N = 163 840, respectively. Galaxies 
were evolved up to i = 12 (3 Gyrs). Results of some 
of the runs are summarized in Table HI where columns 
are: the model label (1), the number of particles (2), 
the SF strength a (3), and SF length scale A (4). As a 
result of simulations the following control parameters are 
displayed: the relative change of components of the disc 
velocity dispersions, measured at time 0.5 and 3 Gyrs (5- 
7) [3^ , the disc angular momentum loss (8), the Toomre's 
Q parameter (9), the Toomre's X parameter (10). The 
expressions for the last two parameters can be found in 
Ref. [13]. 

Table H] shows, as in previous results ^36^, that exper- 
iments made with N = 163 840 is less collisional than 
with smaller number of particles (compare columns 5-7) . 
Runs series SFB are computed with SF strength a = 0.1 
and series SFC with a = 0.3 and we observed that the 
heating of the disc is also higher for a higher SF strength. 

In Fig. 1 we show the time evolution of the amplitude 
of the second harmonic, IA2I, which tells us about the 
appearance of the bar at approximately 1.5 Gyr. The 
bar is stronger for a ~ 0.3 than a = 0.1, meanwhile in 
run SFBOO (Newtonian) a bar appears only at t w 4 Gyrs. 
Also, the disc in presence of SF heats stronger than in 
Newtonian case. This is due to a bar that appears in all 
simulations with SF. 

As it can be seen, a stronger SF produce stronger bars 
for intermediate A, reaching \ A2 \ « 0.3 for models SFC02 
and SFC03. We think that the enhanced heating and 
transfer of the disc angular momentum in runs SFB02 
and SFC03 is due to a stronger and larger bar which in 
turn depends on SF parameters. The bar angular veloc- 
ities neither depend on A nor a and have initially the 
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TABLE I: Numerical parameters of galaxy evolution runs. 



(1) 


(2) 


(3) (4) 


(5) 


(6) 


{') 


(8) 


(9) (10) 


Model 


N 


a 


A 


7r 






-r-^ X 100 


Q 


X2 


SFAOO 


40 960 


0.0 




0.724 0.734 


1.037 


4.1 


2.5 


2.2 


SFAOl 




0.3 


1.0 


1.096 


0.933 


1.404 


6.1 


3.0 


2.0 


SFA02 




0.3 


0.4 


1.051 


0.936 


1.262 


4.5 


2.8 


2.0 


SFA03 




0.3 


0.2 


2.049 


1.576 


1.608 


8.4 


3.0 


1.5 


SFA04 




0.3 


0.1 


1.673 


1.345 


1.114 


7.7 


2.9 


1.0 


SFBOO 


163 840 


0.0 




0.579 


0.447 


0.494 


1.1 


2.0 


3.0 


SFBOl 




0.1 


1.0 0.898 


0.687 


0.636 


2.7 


2.5 


3.0 


SFB02 




0.1 


0.4 


1.176 


0.798 


0.662 


3.3 


2.4 


2.6 


SFB03 




0.1 


0.2 


1.058 


0.805 


0.617 


2.9 


2.3 


2.5 


SFB04 




0.1 


0.1 


0.916 


0.749 


0.624 


2.3 


2.2 


2.7 


SFCOl 


■>■) 


0.3 


1.0 


0.889 


0.653 


0.593 


4.4 


2.5 


2.5 


SFC02 


•>1 


0.3 


0.4 


1.023 


0.785 


0.572 


5.6 


2.6 


2.5 


SFC03 


11 


0.3 


0.2 


1.779 


1.325 


0.764 


8.8 


3.0 


1.7 


SFC04 


11 


0.3 


0.1 


1.279 


1.015 


0.613 


7.7 


2.6 


1.2 



value rjp « 6, which after 1 Gyrs decrease to ftp ~ 5.7 
in code units. 

From Fig. [T] it is seen that SF with a = 0.3 and 
A = 0.2 produces the strongest bar. This is probably 
due to some resonance of disc particle orbits with the se- 
lected A, which is roughly equal to the bar's length. The 
same result can be seen in Fig. [2J where the final density 
contours of discs are plotted. 



IV. INTERACTING GALAXIES SIMULATIONS 

In this section we describe SD-simulations of collisions 
of two identical galaxies, each of which are the same as 
the used in study of isolated galaxies in the preceding 
section. In a recent paper [36j we have shown how the 
numerical parameters influence the bar properties. In 
accordance with this study we choose for the total num- 
ber of particles TV = 163 840, the softening parameter 
e = 0.008 and the time-step At = 1/128. This choice of 
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FIG. 1: Evolution of 1^2! for models SFB00-SFB04 (top 
panel) and for models SFC01-SFC04 (bottom panel). 
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FIG. 2: Projected disc density contours at t = 3Gyrs for 
models SFBOO, SFC01-SFC04. Thick lines indicate zero level 
density contours in logarithmic scale. 
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TABLE II: Parameters of collisions with fixed q — 0.1 



Run 


A 


P 


Wiggle ? Wiggle ? 








Disc 1 


Disc 2 


SFLOfPO 


0.1 





no 


no 


SFLOfPf 


- 


0.4 


yes 


no 


SFL01P2 


- 


0.8 


yes 


yes 


SFL02P0 


0.2 





no 


no 


SFL02Pf 


- 


0.4 


yes 


no 


SFL02P2 


- 


0.8 


yes 


yes 


SFL04P0 


0.4 





no 


no 


SFL04Pf 


- 


0.4 


yes 


no 


SFL04P2 




0.8 


yes 


yes 


SFLIOPO 


1.0 





no 


no 


SFLIOPI 




0.4 


yes 


no 


SFL10P2 




0.8 


no 


yes 


SFNPO 


oo 





no 


no 


SFNPf 




0.4 


yes 


yes 


SFNP2 




0.8 


no 


yes 



parameters prevents from an early bar formation, hence, 
it permits to study the tidal effects on the bar formation. 

We study the effect of the SF on tidal bar formation 
and its properties, such as its amplitude and rotational 
velocity. We compare the bar amplitude and its pattern 
speed for fixed a, varying A and the impact parameter 
(P). 

In order to maintain the same impact velocity and peri- 
centric separation we have studied head-on and off-axis 
impacts of galaxies launched with the initial velocities 
''^ = I^^kI = 200 km/s and the impact parameter p, whose 
values are listed in Table [TTl The galaxies were relaxed 
up to t = 0.25 Gyrs before placed on the orbits with the 
initial separation i? = 64 kpc. The whole collision is fol- 
lowed up to i = 4 Gyrs. We consider prograde-retrograde 
and planar collisions which allow us to investigate two 
possible directions of rotation and to check whether the 
bars emerge in retrograde discs during the violent colli- 
sion. The first galaxy is retrograde, moves to the left and 
for off-axis collisions is placed above, whereas the second 
galaxy is prograde, moves to the right and is located be- 
low the first one. 

The performed collision simulations are summarized 
in Table [TTl where we have varied A for a fixed value 
of the Yukawa strength a = 0.1. We use the follow- 
ing labeling in model names: SF - for the scalar field, 
L(xx) - for the lambda value multiplied by ten, and 
P(y) - for the pericentric parameter expressed in num- 
ber of disc radii. For the Newtonian simulations the 
SFN label is used. In all runs the total energy and 
the total angular momentum conserve better than 1 
%. Movies of some collision simulations are available at 
www. astro . inin.mx/ruslan/stt . 

For all numerical experiments, we have plotted the evo- 
lution of the amplitude of the second harmonic, |^2|, 



which indicates the presence of a bar and corresponding 
pattern velocity, O. We first consider head-on collisions. 
The graphics of | and shown in Figs.[2]and[31 respec- 
tively, are all similar and comparable to the Newtonian 
case, except for small oscillations at the end of the run. 
These oscillations increase with increasing A and are also 
present in plots of bar pattern speed. 

Next, we discuss simulations with an impact param- 
eter equal to the disc's radius, p = 16 kpc. A striking 
difference between the run SFNPl and runs SFLOIPI- 
SFLIOPI is that for SF models the bars in both discs 
have roughly the same amplitude, independently of A, 
whereas in Newtonian case their amplitudes differ by 
roughly twice, see Fig. \E\ . The fact that the retrograde 
discs form bars indicate that the discs in presence of SF 
are unstable and a short and strong enough perturbation 
is sufficient to produce a bar. As in Newtonian case, the 
retrograde bars are slightly faster than the prograde ones, 
indicating their similarity with isolated bars. In general, 
the pattern velocities of the prograde and retrograde bars 
for models with SF are smaller than for the Newtonian 
model, see FiglS] . 

Concerning encounters with p = 32 kpc, the curves in 
Figs. [7| and [8] are similar for each case and show no much 
difference. The only remarkable feature is the higher 
peaks in amplitudes of the prograde discs with SF in 
comparison with the Newtonian case. 

In general, the presence of SF with positive a reduces 
the gravitation force on scales r > A. In Ref. [s^ we 
have seen that larger values of e for the Plummer soft- 
ening produce bars earlier. Thus, it is possible that the 
reduction of gravity either due to small e or large A is re- 
sponsible for early formation of long bars. With negative 
a we expect to have a contrary result, i.e., smaller bars 
that form later. 

In addition, we have looked for morphological differ- 
ences in spiral arms, such as the wiggle effect, which is 
the periodic change of the direction of short arms at the 
end of the bar from leading to trailing [s^ . We visually 
search for the presence of the wiggle in both discs, and 
the results are listed in Table [III 



V. CONCLUSIONS 

We have used the Newtonian limit of general STT 
that are compatible with local observations by the ap- 
propriate definition of the background field constant, i.e., 
< (f) >= G'Ar/(l -f a). Then, from large-scale experiments 
we set a range of possible variations of the parameters of 
the modified gravitational theory, parametrized by (A, a). 
The general gravitational effect is that the interaction 
with the SF becomes weaker {a > 0) by a factor 1/(1 + a) 
for r > A in comparison with the Newtonian case. Using 
the resulting modified dynamics, we have studied isolated 
spirals and the collision of two equal spiral galaxies. From 
our simulations with different A, we have found that the 
inclusion of the SF changes the dynamical properties of 
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FIG. 3: Evolution of jAal for models SFNPO-SFLIOPO. Dot- 
ted lines correspond to the first galaxy (retrograde orbit) 
whereas solid lines correspond to the second galaxy (prograde 
orbit). 



FIG. 4: Evolution of Q. for models SFNPO-SFLIOPO. The 
correspondence of curves is the same as in Fig. O 



galaxies such as the bar morphology and pattern velocity. 

From the performed simulations of isolated galaxy 
models vifith different A, we can see that the addition of 
a non-minimally coupled SF slightly modifies the equi- 
librium of Newtonian model, acting as a perturbation, 
and diminishes the total potential energy, since the ef- 
fective gravitational constant diminishes. This effect 
destabilizes the disc to form a bar in all models with 
iV = 163 840. For a = 0.1, the SF interaction scale 
A = 16 kpc produces a strongest bar, while for a = 0.3, 
a strong bar forms for A = 8 kpc. Also, for these scales 
bars appear earlier. This suggests that there exists some 
kind of resonance between stellar orbits and SF interac- 
tion scale. The results we have found with positive val- 
ues of a imply that most of the spiral galaxies should be 



barred, but this does not exactly correspond with the ob- 
servational fact that around 70% of isolated galaxies are 
barred [45|. However, the bar found might be the result of 
the algorithm to construct the initial models. Therefore, 
a next step is to construct a stable self-consistent model 
in accordance with the modified gravity. On the other 
hand, one could also study the eflFects for negative values 
of a, where the force augments for distances bigger than 
A. A wide range of parameters should be investigated 
and higher resolution have to be used in simulations in 
order to make predictions for particular models. 

In the study of isolated galaxies was shown that the 
presence of the SF destabilizes the disc of isolated galax- 
ies and favors the bar formation. For collisions of two 
galaxies we observe the same trend. In the off-axis colli- 
sions with the impact parameter equal to the disc radius. 
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time, Gyrs 



FIG. 5: Evolution of |^2| for models SFNPl-SFLlOPl. The 
correspondence of curves is the same as in Fig. |31 



the bars in both prograde and retrograde discs have the 
same amplitude, independently of A. However, the wig- 
gle does not appear in the second disc, as shown in Table 
im All these properties depend on the pair (A, a), which, 
on the other hand, can be constrained from observations 
that eventually will discriminate among the different val- 
ues of the parameters of the theory. 

The results presented are only preliminary, and we de- 
scribe the overall differences without giving a full inter- 
pretation. A broad range of parameters should be inves- 
tigated and higher resolution have to be used in simu- 
lations in order to make comparisons with the observed 
interacting galaxies. 
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the potential oscillates around the background field 

(0)= const., (A.l) 

and we expand the field quantities around (0) and the 
Minkowski metric (g^jy = r]^i, + h^i,) to first order, being 
{(p) « (f> and hfj^i, « r]^v We assume that the field 
is quasi-stationary, such that all time derivatives can be 
ignored: = — V^0. Additionally, one has that p(? » 
p, and T = Too = p. 

With these assumptions we can expand all terms of 
field equations (013]) and use Taylor series for the un- 
known functions uj{4>) and V{(j)) in terms of the small 
quantity (f>. Let us first consider Eq. The terms on 
the l.h.s., after inserting the perturbed metric and as- 
suming the linearized harmonic gauge, can be written 
as: 



(A.2) 



W,. The first term on the r.h.s. is the usual 



where h 

general relativity term. The second term can be written 
as follows: 



(A.3) 



The third and fourth terms vanish because these are sec- 
ond order terms in the expansion: 



da 



5u 



The last two terms can be written as 



and 



0- 
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(A.5) 
(A.6) 



Substituting the above expansions into Eq. ^ we may 
write: 
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APPENDIX: THE NEWTONIAN LIMIT OF STT 

We start with the field equations are given bye Eqs. 
and ([3]) and compute the Newtonian limit of the scalar- 
tensor theory to first order. Accordingly, we assume that 



V 



1 



(0) {^y 



(A.7) 



If we multiply the Eq. (jA.7p by Ty^'' and take the trace, 
the l.h.s. can be written as 



1 



and the r.h.s. as 



V 



2" 
(A.8] 



{4>) (0)' 



87r/9 - -(4 + /i) - V^(t> + (4 + h)V'^(, 



(A.9) 
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In the last equation we used the relation 



(A.IO) 



valid to first order. Now, let us rewrite the Eq. ()A.2p as 
follows: 



1 



1 



in which we insert IA.9I instead of the term — and 
using Eq. (jA.7|) we obtain 



V 



Sirp - -(4 + h) + hV^4> + 



{4>) {c^y 



V 



(A.12) 



Let p = V = Q. Then, the above equation becomes to 
first order 



2 "00, A 



Vh 1. 



V V V,. X o- 

47rpH /loo H V"^' 

2 2 " 4 2 







v 









(A.13) 



where $ is modified Newtonian potential. Here we used 
the fact that 

\ = ^ (^oo)'"o - \ = {Koy\^ = 



1 



Finally, expanding the potential 

dV 

V = vm) + -r 



-V^/ioo. (A.14) 



(A.15) 



(0) 



and discarding terms higher than first order we reach the 
the expression 



1 



(A.16) 



Accordingly, particles in our galactic models experience 
a force given by F = — V$jv- Now, consider Eq. ([3]). 
The terms that depend on can be expanded to the first 
order: 



1 



1 



3 + 2lj 3 + 2lj 



1 



d 



d(j) \3 + 2uj 



(0) 



Let us make the following notation: 

dV 



1 



3 + 2u} 



i'p)' 



V, 



(0) 



d 



1 



3 + 2w 



(A.17) 
(A.18) 



(A.19) 



(A.20) 

Substituting Eqs. ((XT5| . (|A.17IIA.20p into Eq. ^ and 
keeping only the terms up to the first order, we obtain: 



{ml 



{<p) 



2Vi 



(0) 



87Tp + 2Vi^^^~{m 



(0) 



(A.21) 



The constant in the second term of the l.h.s. represents 
an -squared- effective mass term of the theory that we 
will denote as m^. Last two terms in the r.h.s. represent 
a cosmological constant that we will set to zero, since 
the mean density of the galaxy is much greater than this 
term. Accordingly, we finally have: 



(A.22) 



This equation together with Eq. (|A.16p represent the 
Newtonian limit of general STT that can be expanded 
around the background quantities (0) and 77^^. 



